Structural insights into TRPV2 activation by small molecules

Transient receptor potential vanilloid 2 (TRPV2) is involved in many critical physiological and pathophysiological processes, making it a promising drug target. Here we present cryo-electron microscopy (cryo-EM) structures of rat TRPV2 in lipid nanodiscs activated by 2-aminoethoxydiphenyl borate (2-APB) and propose a TRPV2-specific 2-ABP binding site at the interface of S5 of one monomer and the S4-S5 linker of the adjacent monomer. In silico docking and electrophysiological studies confirm the key role of His521 and Arg539 in 2-APB activation of TRPV2. Additionally, electrophysiological experiments show that the combination of 2-APB and cannabidiol has a synergetic effect on TRPV2 activation, and cryo-EM structures demonstrate that both drugs were able to bind simultaneously. Together, our cryo-EM structures represent multiple functional states of the channel, providing a native picture of TRPV2 activation by small molecules and a structural framework for the development of TRPV2-specific activators.

T ransient receptor potential vanilloid 2 (TRPV2) is a Ca 2+permeable non-selective cation channel belonging to the vanilloid subfamily of the transient receptor potential (TRP) channels 1,2 . TRPV2 is expressed in most organs of the human body, including the heart, white blood cells, the pancreas, the central and peripheral nervous systems, the bladder, the prostate, and the placenta 3 . TRPV2 has been shown to play a significant role in maintaining physiological cardiomyocyte function, with a potential therapeutic use of TRPV2 blockade in cardiomyopathy [4][5][6] . TRPV2 plays a role in phagocytosis in macrophages 7,8 , placental development 9 , T-cell activation 10,11 , and insulin secretion in pancreatic β-cells 12,13 . Beyond its routine functions in healthy cells, TRPV2 has been shown to play a significant role in the progression, drug resistance and metastasis of different forms of cancer 14 . Thus, TRPV2's role in both healthy and malignant tissues makes it a promising target for drug discovery and the development of novel therapeutics in cardiovascular diseases, cancer, and additional pathophysiological conditions.
A key limiting factor in understanding precise TRPV2 molecular function is the lack of specific modulators of the channel. While endogenous modulation of TRPV2 remains largely unexplored, there are several known non-specific exogenous activators of TRPV2, including 2-aminoethoxydiphenyl borate (2-APB), cannabidiol (CBD) and probenecid, and exogenous inhibitors, including ruthenium red, trivalent cations, tranilast, and piperlongumine 2,15 . To aid the rational development of modulators with increased efficacy and TRPV2 specificity, more structural information is needed to determine binding sites for the known ligands and track channel gating.
Early structural studies of TRPV2 in the absence of TRPV2 modulators hypothesized mechanisms of channel gating [16][17][18][19][20] , but due to modifications in the protein sequence and the use of detergents those structures did not provide a comprehensive picture of wildtype TRPV2 gating in a native lipid environment. Recently, our lab reported CBD-bound TRPV2 structures in lipid nanodiscs revealing CBD binding to a therein reported ligand binding site 21 . Nevertheless, the CBD-bound TRPV2 structures did not reveal a mechanism for TRPV2 opening by small molecules 21 . Thus, it is essential to study how other known TRPV2 exogenous activators interact and gate the channel in order to get a full picture of the gating cycle and discover potential ligand-binding sites.
To address the need for structural information of TRPV2 ligand-binding and gating, here we use cryo-electron microscopy (cryo-EM) to determine the structures of wildtype rat TRPV2 in the presence of 2-APB in lipid nanodiscs, which reveals a binding pocket for 2-APB between the S5 helix and the S4-S5 linker of two adjacent monomers. In silico docking and electrophysiological data validate the proposed binding site for 2-APB and confirm a key role for His521 and Arg539 in the binding pocket. Additionally, electrophysiological data reveal that CBD acts as a potentiator for 2-APB activation of TRPV2 and cryo-EM structures confirm that both drugs can bind simultaneously at their unique binding pockets. Cryo-EM also allow us to resolve two distinct conformations of TRPV2 bound to 2-APB: an activated state and an inactivated state, leading us to propose a wedge mechanism for TRPV2 modulation by 2-APB. While 2-APB does not bind in the vanilloid pocket, it induces conformational changes in TRPV2 that lead to the expulsion of the vanilloid lipid as seen in other TRPV channels upon activation. Together, our results put TRPV2 in line with activation mechanisms established for other TRPV channels and lay the foundation for future development of TRPV2specific modulators.

Results
Activated and inactivated structures of 2-APB-bound TRPV2. 2-APB is a well-established activator of the TRPV1-3 channels 22 . Recent structural studies 23-25 revealed a consensus 2-APB binding site for TRPV3 positioned between the Pre-S1 helix and the TRP helix, which agrees with previous functional data 26 . Since the critical histidine residue of the TRPV3 2-APB binding site is not conserved in TRPV2 26 ( Supplementary Fig. 1A), we hypothesized that TRPV2 interacts with 2-APB through a different binding site that could be used to develop TRPV2-specific therapeutics and unravel the TRPV2 gating mechanism.
In order to obtain an activated structure of the TRPV2 channel and identify the TRPV2 2-APB binding site, we collected a large cryo-EM dataset of rat TRPV2 reconstituted in lipid nanodiscs and incubated with 1 mM 2-APB for 5 min before preparing the grids. The dataset showed extensive heterogeneity and was analyzed using the 3D variability tool 27 from cryoSPARC 28 ( Supplementary Fig. 2). This generated continuous families of structures reflecting the major movements that appear in the sample and assigned each particle to a specific position within each family. The family of states for the C4 mode of movement for 2-APB-bound TRPV2 showed substantial conformational changes and partial opening of the channel (Supplementary Movie 1). The particles corresponding to the most open and the most closed states were isolated and refined independently, providing two distinct structures with C4 symmetry which we named TRPV2 2-APB-bound "activated" state (TRPV2 2APB_AC ) at 4.15 Å resolution and TRPV2 2-APB-bound "inactivated" state (TRPV2 2APB_IAC ) at 3.5 Å resolution ( Supplementary Fig. 3, Table 1).
We compared TRPV2 2APB_AC and TRPV2 2APB_IAC with each other and with the previously identified fully closed apo state of the channel in nanodiscs (TRPV2 Apo1 , PDB 6U84) 21 and all three structures were remarkably different at their pore profiles and at the cytoplasmic regulatory switch 23 (Fig. 1A-D, G-I). The pore of TRPV2 2APB_AC was partially open (Fig. 1B) with a pore radius of 2.0 Å at Gly606 in the selectivity filter and 1.9 Å at the Met645 in the lower gate (Fig. 1D). Thus, the diameter of the TRPV2 2APB_AC pore could allow a partially dehydrated Ca 2+ ion to go through. The S6 helix of TRPV2 2APB_AC remained entirely α-helical similar to TRPV2 Apo1 (Fig. 1A, B), but with a significant conformational change at the bottom of S6 where the sidechain of His651 was rotated out of the pore by 180°. In this new position, His651 stabilized TRPV2 2APB_AC in the partially open conformation through the newly formed cation-π interaction between His651 and Lys531 from the S4-S5 linker (Fig. 1B, E). In contrast, in TRPV2 Apo1 His651 tightly seals the lower gate (Fig. 1A, E). At the cytoplasmic regulatory switch, TRPV2 2AP-B_AC was similar to TRPV2Apo1 (Fig. 1G) and characterized by the presence of a stable C-terminal helix (residues 710-720, Fig. 1H) and N-terminal loop (residues 30-45, Fig. 1H), suggesting that we captured TRPV2 in an activated state that is an intermediate or partially open.
The TRPV2 2APB_IAC state was closed at both lower and upper gates (Fig. 1C, D). The narrowest constrictions in the pore were at Gly606 in the selectivity filter and at the Met645 in the lower gate, where the pore radius was 0.6 Å and 0.9 Å, respectively (Fig. 1C, D). The S6 helix of TRPV2 2APB_IAC was entirely α-helical (Fig. 1C) and the His651 at the bottom of S6 was still rotated out of the pore, forming a cation-π interaction with Lys531 ( Fig. 1F). At the cytoplasmic regulatory switch, TRPV2 2APB_IAC had a C-terminal loop (residues 710-726, Fig. 1I) which displaced the N-terminal loop. Therefore, we interpreted TRPV2 2APB_IAC as an inactivated state of TRPV2.
Identification and validation of the 2-APB binding site. The TRPV2 2APB_IAC structure revealed a putative 2-APB binding site between S5 of one TRPV2 monomer and the S4-S5 linker of the adjacent ( Fig. 2A, B, Supplementary Fig. 4). The proposed TRPV2 2-APB binding site is located at the N-terminal end of the S4-S5 linker, where the ligand coordinates with Arg539, Thr522, His521, and Tyr525 ( Fig. 2A, B). In the rat TRPV2 Apo1 21 , Arg539, His521, and Tyr525 form a network of cation-π interactions, which stabilize the interaction between the S4-S5 linker and S5 from an adjacent monomer (Fig. 2C). 2-APB binding disrupts this network and may act as a wedge to force apart two adjacent TRPV2 monomers at the interface between the S4-S5 linker and the bottom of S5 ( Fig. 2A, B, Supplementary Fig. 4). To validate the position and orientation of the 2-APB molecule in the binding site we used RosettaLigand docking 29,30 . The top-scored docking model was in the same location and similar orientation as in the initial fitted model with the reported all-atom root-mean-square deviation (RMSD) of 2.106 Å between the two positions of the drug (Fig. 2D, Supplementary Fig. 5).
To further validate this binding site for 2-APB, we performed whole-cell patch clamp recordings on HEK293T cells transiently transfected with either the wildtype rat TRPV2 or with TRPV2 mutants carrying single-point mutations around the proposed binding site: His521Ala, His521Asn, Thr522Ala, Tyr525Ala, Arg535Lys, and Arg539Lys. As histidine and arginine residues were critical for the TRPV3 consensus 2-APB binding site, we followed the previously established substitutions for these residues 26 . Increasing concentrations of 2-APB, from 30 up to 5000 µM, were applied to the cells in order to obtain dose-response curves for 2-APB-evoked activation for each construct (Fig. 2E-J, Supplementary Fig. 6). The calculated EC 50 values are summarized in Table 2. When compared to the wildtype TRPV2 (Fig. 2E, Table 2), the 2-APB sensitivity in the His521Ala, His521Asn, and Arg539Lys was reduced ( Fig. 2F-H, Table 2) and the double mutant His521Ala/Arg539Lys showed a fivefold decrease in 2-APB sensitivity (Fig. 2I, Table 2). This reduction is comparable to the loss of 2-APB sensitivity previously reported for TRPV3 upon mutation of His426 or Arg696 26 . All mutants involving His521 and Arg539 responded normally to high heat ( Supplementary Fig. 6), demonstrating that the mutated channels were functional and that the effect of these mutations was specific to 2-APB. Interestingly, the Thr522Ala and Arg535Lys mutations caused a mild increase in 2-APB sensitivity ( Supplementary Fig. 6, Table 2). Finally, Tyr525Ala seemed to have a global loss-of-function phenotype as we could not evoke any membrane currents with either 2-APB or heat. These data strongly support a potential interaction of 2-APB and rat TRPV2 via the proposed binding site and suggest a key role of His521 and Arg539 in the binding.
The 2-APB binding site determined here for rat TRPV2 is unique to TRPV2 channels. Interestingly, it features the same key amino acids as the consensus 2-APB binding site in mouse TRPV3-histidine and arginine 26 . His521 is highly conserved in mammalian TRPV2, including human TRPV2 ( Supplementary  Fig. 1B), supporting the recent electrophysiology data on human TRPV2 activation by 2-APB 31 which differs from earlier calcium flux studies 32,33 . Notably His521 is not conserved in other human TRPV channels ( Supplementary Fig. 1A). This mirrors the uniqueness of the 2-APB binding site in TRPV3, which features a conserved His426 for mammalian TRPV3 26 ( Supplementary  Fig. 1A). Together these results are in full agreement with Simultaneous binding of 2-APB and CBD to TRPV2. Given the proximity of our herein proposed 2-APB binding site to our previously determined CBD binding site 21 , we decided to investigate whether these two activators could bind simultaneously and affect channel activation. Indeed, when applied together, 2-APB and CBD have a synergetic effect on TRPV2 activation as shown by whole-cell patch recordings on HEK293T cells transiently transfected with the wildtype rat TRPV2 ( Fig. 3A-C). Application of 2-APB alone evoked small currents with a density of 0.98 ± 0.1 pA/pF, 2.8 ± 0.3 pA/pF, and 24.9 ± 3.7 pA/pF at 10, 30, and 100 μM respectively (Mean ± SEM, N = 7, Fig. 3A). At the same time, in combination with 10 μM CBD 3 μM 2-APB was already enough to elicit a current with a density of 8.9 ± 6.6 pA/pF (Fig. 3B); the current density upon co-application of 10 μM CBD with 10, 30, and 100 μM 2-APB was equal to 51.4 ± 21.3 pA/pF, 189.4 ± 21.1 pA/pF, 260.9 ± 29.5 pA/pF respectively (Mean ± SEM, Fig. 3B). Thus, a combination of two drugs produced ten times greater current density than the current density induced by either of the agonists alone ( Fig. 3B, C). Therefore, we next collected a cryo-EM dataset of rat TRPV2 in lipid nanodiscs incubated with 1 mM 2-APB and 30 μM CBD for 30 min before preparing the grids. The double-drug dataset had similar heterogeneity to that observed in the 2-APB dataset and was again analyzed using the 3D variability tool 27 from cryoSPARC 28 (Supplementary Fig. 7). The dataset yielded two states with C4 symmetry: one identified as activated based on the conformation of the cytoplasmic regulatory switch 23 with a stable C-terminal helix and N-terminal loop at 3.7 Å resolution (TRPV2 2APB_CBD_AC ) and the other identified as inactivated with a C-terminal loop at 3.8 Å resolution (TRPV2 2APB_CBD_IAC ) ( Supplementary Fig. 8, Table 1  . The radius of a dehydrated calcium ion is marked by a dotted gray line at 1.1 Å. E, F His651 is inserted into the pore in the apo state (TRPV2 Apo1 , wheat), but rotates out of the central channel of the pore on channel activation (TRPV2 2APB_AC , light blue) to form a cation-π bond with Lys531. During channel inactivation (TRPV2 2APB_IAC , light cyan), His651 remains stably rotated out of the pore. TRPV2 Apo1 (G) and TRPV2 2APB_AC (H) feature the active conformation of the cytoplasmic regulatory switch, with the C-terminus forming a helix (green) inside the ankyrin repeat domain (ARD) skirt and the N-terminus (yellow) forming a loop at the multimerization interface between the β-sheet region and an ARD. TRPV2 2APB_IAC (I) features the inactive conformation of the cytoplasmic regulatory switch, with the C-terminus (red) forming a loop that wraps around the β-sheet region.
site 21 (Fig. 3B, C, Supplementary Fig. 4). The drug densities were well-resolved in the TRPV2 2APB_CBD_AC structure, where we were able to build both 2-APB and CBD in orientations similar to the poses found in the single drug structures (Fig. 3B, C) 21 . Therefore, the double-drug data independently supported the proposed 2-APB binding site. Our previously published CBD-bound TRPV2 CBD1 structure 21 , which has an "activated" conformation of the cytoplasmic regulatory   Fig. 9; 0.751 Å Cα RMSD between TRPV2 2APB_IAC and TRPV2 CBD2 ; 0.761 Å Cα RMSD between TRPV2 2APB_CBD_IAC and TRPV2 CBD2 , PDB 6U88) and only differed locally at the drug binding sites. This led us to the conclusion that transition to the activated conformation in both datasets is driven by 2-APB binding rather than CBD binding, with 2-APB acting as a wedge to initiate conformational changes in the channel.
We next analyzed the pore profiles of TRPV2 2APB_CBD_AC and TRPV2 2APB_CBD_IAC (Fig. 3F-H). Both double-drug structures remained closed at the selectivity filter; the lower gate of the TRPV2 2APB_CBD_AC was partially open with a pore radius of 1.2 Å at Met645 compared to 0.3 Å in the TRPV2 Apo1 and 0.5 Å in the TRPV2 CBD1 (Fig. 3H, Supplementary Table 1). The TRPV2 2APB_CBD_IAC structure was closed at Met645 with a pore radius of 0.7 Å, similar to TRPV2 CBD2 at 0.8 Å (Supplementary  Table 1). However, these radii are likely insufficient for permeation of calcium ions through the hydrophobic pore of TRPV2. In both double-drug structures, the His651 sidechain was rotated out of the pore as in the 2-APB-bound structures.
The mechanisms of TRPV2 activation and inactivation. To understand the gating cycle upon 2-APB binding, we next compared the conformational changes observed between TRPV2 Apo1 and TRPV2 2APB_AC ( , aligned with respect to the pore-forming S6 helices (residues 618-648). These assessments revealed that TRPV2 gating could be described as the opening and closing of an iris aperture. To transition from the TRPV2 Apo1 closed state to the TRPV2 2APB_AC open state, the entire assembly except the pore domain (S5-PH-S6) rotated clockwise around the central pore axis by 7°(if viewed from the cytoplasmic side of the protein) (Fig. 4A, B). This large movement of the entire channel originated from a 0.8 Å shift of the S4-S5 linker away from the central pore axis (Fig. 4B). In the pore domain, S5 and the pore helix shifted upwards by 1.4 Å (Fig. 4B), leading to the opening of the selectivity filter together with rotation of Met607 out of the pore (Fig. 1A, B). The expansion of the lower gate was achieved through the widening of the entire transmembrane domain by 0.8-1.0 Å (Fig. 4C), which led to a rotation of His651 at the bottom of S6 out of the pore (Fig. 1A, B). Finally, the ankyrin repeat domains (ARDs) show a slight downward movement of the channel by 1.1 Å (Fig. 4C).
On the contrary, to transition from the activated to the inactivated state, the entire assembly except the pore domain (S5-PH-S6) rotated counterclockwise around the central pore axis by 7°( if viewed from the cytoplasmic side of the protein) (Fig. 4D, E). The S4-S5 linker and the TRP helix rotated and simultaneously shifted by 1.5 Å towards the central pore axis (Fig. 4E) as part of the narrowing of the entire transmembrane domain of the channel (Fig. 4F) which closed the lower gate of the pore (Figs. 4F, 1C, 1F). The selectivity filter closed through a downward shift of S5 and the pore helix by 1.4 Å (Fig. 4F). Finally, in the inactivated state there was a significant 4.4 Å downward movement of the ARDs (Fig. 4F), caused by the combination of a downward shift of 2.3 Å and a 5°p ivot of each individual ARD (Fig. 4F). The ARDs pivoting coincided with the conformational changes in the cytoplasmic regulatory switch (Fig. 1H, I). In order to come back to the TRPV2 Apo1 closed state from the 2-APB-bound TRPV2 2APB_IAC inactivated state, the drug needs to leave the binding pocket allowing the structure to relax the entire transmembrane domain and ARDs to their original positions (Fig. 4G-I).
The double-drug structures displayed similar movements throughout the gating cycle (Fig. 3G, H, Supplementary Fig. 9). Comparing to TRPV2 Apo1 , the TRPV2 2APB_CBD_AC state displayed the same 7°rotation around the central pore axis as TRPV2 2APB_AC and a 1.7 Å shift of the S4-S5 linker away from the central pore (Fig. 3G). While the shift of the S4-S5 linker was larger in TRPV2 2APB_CBD_AC than in TRPV2 2APB_AC , this is likely the result of local accommodation of CBD as the overall expansion of the transmembrane domain did not increase. The other significant difference between the two structures was observed at the top of the pore, where in the double-drug structure the pore helix did not move up as much as in the 2-APB-bound one (Fig. 3G) and therefore the selectivity filter of the double-drug structure remained closed (Fig. 3D, F). This is primarily because of the presence of CBD, which stabilizes S5 and does not allow for its upward movement and thus for propagation of conformational changes upon 2-APB binding all the way up to the pore helix. The inactivated structures TRPV2 2APB_IAC and TRPV2 2APB_CBD_IAC were mostly identical and only differed at the drug-binding sites with small local adjustments (Fig. 3H). Therefore, we concluded that 2-APB-bound and double-drug structures represent the same functional states of the channel.
The double-drug structures were also complimentary to the 2-APB structures when it comes to analysis of the lipid density in the vanilloid pocket as part of the TRPV2 activation mechanism. In the closed TRPV2 Apo1 state, as well as in TRPV2 Apo2 21 and TRPV2 CBD1 21 , the vanilloid pocket contained well-resolved lipid tails with no definitive head group (Fig. 5A, Supplementary  Fig. 10). We compared these lipid tail densities with the corresponding non-protein densities in the vanilloid pockets of the best 2-APB-bound structures: activated TRPV2 2APB_CBD_AC and inactivated TRPV2 2APB_IAC (Fig. 5). In TRPV2 2APB_CBD_AC , the lipid underwent significant conformational re-arrangement with the lipid resting along the S4-S5 linker rather than extending upwards along S4 (Fig. 5B, Supplementary Fig. 10). On the other hand, TRPV2 2APB_IAC had minimal lipid density, similar to our previously solved inactivated CBD-bound structure TRPV2 CBD2 21 , and the cryo-EM map was of sufficient quality to resolve the lipid if it were present (Fig. 5C, Supplementary  Fig. 10). The other two 2-APB-bound states, TRPV2 2APB_AC and TRPV2 2APB_CBD_IAC , were consistent with the structures described above (Supplementary Fig. 10). Therefore, we concluded that the vanilloid lipid may play a key role in 2-APB-based TRPV2 channel gating, similar to the previously proposed mechanism for TRPV1 activation by small molecules 34 .

Discussion
In the present study we were able to determine structures of 2-APB-bound TRPV2 embedded in lipid nanodiscs in two distinct states, which provided a unique opportunity to analyze modulation of this channel by 2-APB in a near-native environment. TRPV1-TRPV3 are structurally homologous, share at least 38% sequence identity and can be activated by 2-APB 22 , thus it is logical to assume that they share a common 2-APB interaction site. However, there is compelling evidence that 2-APB binding sites could be unique to each channel. In 2009, the Patapoutian lab used a high-throughput mutagenesis screen of mouse TRPV3 to identify His426 and Arg696 as crucial for effective activation by 2-APB 26 . Notably, His426 is not conserved in TRPV1 or TRPV2 and introducing it at that location did not change the effect of 2-APB on either channel 26 . In 2010, the Vlachova lab used targeted mutation of the rat TRPV1 S4-S5 linker to identify mutations that reduced 2-APB sensitivity without loss of heat or capsaicin sensitivity, including Lys571Glu 35 . They made the analogous mutations in TRPV2 and TRPV3 and did not observe any loss of 2-APB sensitivity 35   previously characterized His426/Arg696 site 25 . The other two binding sites were not observed in the subsequent cryo-EM studies of 2-APB-bound TRPV3 23,24 , nor in the 2-APB-bound TRPV2 structures presented in this paper ( Supplementary  Fig. 11). However, one of these additional 2-APB binding sites was identified as the 2-APB binding site in TRPV6 where 2-APB acts as an inhibitor 36 . Overall, these previous studies show that the interaction of 2-APB with TRPV1-TRPV3 is not based around a single shared site [23][24][25][26]35 .
Our cryo-EM data allows us to propose a single novel site for 2-APB binding to TRPV2 between S5 and the S4-S5 linker of an adjacent monomer (Fig. 2). In silico docking supports the putative 2-APB binding site and the orientation of the drug. Using electrophysiology, we verified the key interacting residues by mutations but were not able to completely abolish TRPV2 sensitivity to 2-APB. Notably, the previous studies to identify the 2-APB binding site in TRPV3 were also not able to completely abolish the response to 2-APB 25,26 . This may indicate that TRPV2 contains more than one 2-APB binding site, as was proposed for TRPV3 by the Sobolevsky lab, but we observe only one potential location for 2-APB in our current cryo-EM maps.
To date, our TRPV2 structures in the presence of CBD 21 or 2-APB in a lipid environment have revealed two unique locations for activator binding at the dynamic S4-S5 linker region, sometimes referred to as the "gearbox" of TRP channels 37 . This data argues for the S4-S5 linker as a rich target area for development of both TRPV2-specific activators and inhibitors. Despite the close proximity of the binding sites, the two drugs can bind to TRPV2 simultaneously (Fig. 3) and function synergistically. The addition of CBD appears to have a stabilizing effect on the activated state as seen from the larger proportion of particles and improved resolution of the activated structure in the double-drug dataset. The proposed mechanism of 2-APB action as a wedge further narrows down the region of interest for the rational design of TRPV2-specific drugs to the N-terminal part of the S4-S5 linker.
Abundant cryo-EM and X-ray crystallography data from TRPV1-TRPV3 has become available in the last decade [15][16][17][18][19][20][21][23][24][25]34,[38][39][40][41][42][43][44][45][46][47][48] , but comparison of those channels and identification of structural features functionally relevant for native channel gating remain difficult as different orthologous proteins with sequence modifications were used and the majority were purified in detergents/amphipols for structural studies. Recently, the structures of full-length TRPV1-TRPV3 have been solved in nanodiscs 19,21,24,43,[46][47][48] closely mimicking the native membrane environment and providing the opportunity for their most accurate comparison. Based on this wealth of information, there are a few key common features observed in the gating cycle that are conserved between TRPV1-TRPV3. First, the cytoplasmic regulatory switch at the C-terminus of the channel must be in a helical conformation for channel activation, while a loop conformation of the C-terminus prevents channel opening 23,46,47 (Supplementary Fig. 12). Second, channel activation includes a 4-8°rotation of the ARD and TRP helix around the central pore [46][47][48] (Supplementary Fig. 12), while channel inactivation coincides with a downwards movement of the ARDs 46,47 . Third, the selectivity filter of TRPV1-TRPV3 is not an activation gate and can flexibly accommodate a range of ions 48,49 .
In the current study, we obtained structures of native rat TRPV2 in lipid membranes in the presence of the activator 2-APB alone or together with CBD, resolving a gating cycle of TRPV2. The global movements upon 2-APB activation of TRPV2 are the same as observed in drug and heat activation of TRPV1 and TRPV3: 2-APB binding wedges apart two monomers of TRPV2, rotating the channel 7°around the central pore, disrupting the vanilloid lipid, and initiating pore opening. In the inactivated states, the channel rotated back by 7°and the cytoplasmic regulatory switch transitioned from a helix to a loop conformation, allowing the ARDs to move downwards and closing the channel, but the drug remained bound and the vanilloid lipid was absent.
While the global channel movements of the activated TRPV2 states align with those of the most open TRPV1 and TRPV3 structures ( Supplementary Fig. 12), we also need to look at the pore to determine whether or not the channel is open. The selectivity filter in the two activated states varies (Figs. 1, 3), presumably due to the inherent flexibility of this region of TRPV2. The pore radius of the channel at the lower gate of the TRPV2 2APB_AC is 1.9 Å, which would allow passage of a dehydrated ion (1.1 Å for calcium), but not for a fully hydrated ion (4.1 Å for calcium). The hydration state of ions passing through TRPV1-TRPV3 is unclear, but several structures of TRPV1 and TRPV3 have a lower gate with a minimum radius greater than 1.1 Å (Supplementary Fig. 13, Supplementary Table 1) and their interpretation falls along a spectrum from "sensitized" or B TRPV2 2APB_CBD_AC has disrupted lipid density in the vanilloid pocket. C TRPV2 2APB_IAC displays minimal density that could be attributed to the lipid tails in the vanilloid pocket. One monomer is colored wheat, medium purple or light cyan, the adjacent monomer is colored gray; 2-APB is colored orange, CBD is colored green, potential lipid densities are colored khaki. Maps contoured at σ = 5.
"intermediate" to fully open. Based on a comparison to the other published TRPV1 47,48 and TRPV3 46 structures, we interpret the 2-APB-activated TRPV2 structures as intermediates on the path to opening. Nevertheless, the pore radius in even the most open of the TRPV1-TRPV3 structures is still too narrow to allow fully hydrated ions to pass 50 .
With the extensive progress in cryo-EM imaging and data processing achieved over the last few years, more work can now be done to characterize native modulators of channel activity and for drug discovery.

Methods
Protein sequence alignment of the TRPV channels. TRPV sequence alignments were performed using Clustal Omega 51 . Graphics of sequence alignments were generated using Aline 52 .
Protein expression and purification. Full-length rat TRPV2 was expressed and isolated as previously published 21 . Briefly, rat TRPV2 cloned into a YepM vector and tagged with 1D4 was expressed in Saccharomyces cerevisiae. The yeast expressing TRPV2 were lysed in 25 mM Tris, pH 8.0, 300 mM Sucrose, and 5 mM EDTA using a M110Y Microfluidizer (Microfluidics), and the membranes isolated by ultracentrifugation at 100,000 x g. These TRPV2-enriched membranes were solubilized in 20 mM Hepes, pH 8.0, 150 mM NaCl, 5% glycerol, 0.087% LMNG, 2 mM TCEP, and 1 mM PMSF for 1 h. This mixture was clarified by centrifugation at 100,000 x g and the insoluble fraction discarded. The solubilized TRPV2 was bound to 1D4 antibody coupled CnBr-activated Sepharose beads, followed by washes with Wash Buffer (20 mM Hepes, pH 8.0, 150 mM NaCl, 2 mM TCEP) supplemented with 0.006% (w/v) DMNG. The protein was eluted with Wash Buffer supplemented with 0.006% (w/v) DMNG and 3 mg/mL 1D4 peptide. The purified TRPV2 was concentrated to a volume under 1 mL and reconstituted into nanodiscs in a 1:1:200 ratio of TRPV2:MSP2N2:soy polar lipids (Avanti). The soy polar lipids were rapidly dried under a nitrogen flow and further dried under vacuum before being resuspended in Wash Buffer containing DMNG in a 1:2.5 ratio (soy polar lipids:DMNG). The assembled nanodisc reconstitution mixture was incubated at 4°C for 30 min before adding Bio-Beads to the mixture. After 1 h, the reconstitution mixture was transferred to a new tube with fresh Bio-Beads and incubated overnight at 4°C. The nanodisc embedded TRPV2 was purified from empty nanodiscs using a Superose 6 column (GE) equilibrated in Wash Buffer for size-exclusion chromatography. The eluted TRPV2 was concentrated to 2 mg/mL to use in vitrification.
MSP2N2 was expressed and purified as previously described 53 . Briefly, MSP2N2 inserted into a pET28a vector (Addgene) was expressed in BL21 (DE3) cells. After harvest, the cells were resuspended in a buffer containing 20 mM Tris-HCL, pH 7.4, 1 mM PMSF, and a complete EDTA-free protease inhibitor cocktail tablet (Roche) and lysed using a M110Y Microfluidizer (Microfluidics). The lysate supernatant was bound to Ni-NTA resin, which was then washed with Wash Buffer (20 mM Tris-HCl, pH 7.4, 100 mM NaCl), then by Wash Buffer supplemented with 1% Triton X-100, then by Wash Buffer supplemented with 50 mM sodium cholate, then by Wash Buffer supplemented with 20 mM imidazole. MSP2N2 was then eluted from the Ni-NTA resin with Wash Buffer supplemented with 300 mM imidazole and exchanged into 50 mM Tris-HCl, pH 7.5, 100 mM NaCl, 5 mM EDTA before being concentrated to~10 mg/mL.
Cryo-EM sample preparation. For 2-APB-bound TRPV2, 1 mM 2-APB was incubated with TRPV2 for 5 min on ice prior to grid preparation. For the 2-APB and CBD-bound sample, 1 mM 1-APB and 30 μM CBD were incubated with TRPV2 for 30 min prior to grid preparation. Immediately before blotting, 3 mM fluorinated Fos-choline 8 (Anatrace) was added to improve particle distribution in vitreous ice. 3 μL of sample was applied to a freshly glow discharged 200 mesh Quantifoil 1.2/1.3 grid and then blotted for 5-8 s at 4°C and 100% humidity before vitrification in liquid ethane in Vitrobot Mark IV (FEI).
Cryo-EM data processing. For the 2-APB dataset, the movies were motion corrected the using the Relion 3.1 54 implementation of MotionCor2, binning the pixels to 1.06 Å/pix. The remaining processing was done using cryoSPARC v3.2. CTF fit of the motion corrected micrographs was estimated using patch CTF. The dataset was split into two sets due to the large size-15,237 micrographs in set 1, 19,579 micrographs in set 2. The micrographs in each set were curated to remove suboptimal micrographs, resulting in 13,597 micrographs in set 1 and 18,768 micrographs in set 2. From these sets, template-based auto-picking yielded an initial set of 1,928,671 and 2,517,442 particles, respectively. These particles were binned by 4 during extraction and subjected to multiple rounds of 2D classification to remove false positives and bad particles, reducing the particle count to 732,794 particles in set 1 and 1,125,082 particles in set 2. Both sets of particles were reextracted binned by 2 and subjected to multiple rounds of heterogeneous refinement with the TRPV2 Apo1 map (EMD-20677) as a good reference class and an unrelated map as a reference for noise classes. For set 1/2, this removed 192,713/ 239,710 particles that did not align to TRPV2 classes and 138,471/258,949 particles that aligned to a TRPV2 class with a completely disrupted pore domain, leaving 401,610/626,423 particles in the classes with the best features. These particles were reextracted to their final box size and refined by Non-uniform refinement with C4 symmetry before symmetry expansion. These particles were then subjected to 3D variability analysis 27 , using a mask that excluded nanodisc density and searching for three modes of movement. The resulting modes were primarily assessed based on the cytoplasmic regulatory switch, as it is a well resolved region with a binary variable. Mode 1 from set 1 and Mode 2 from set 2 represented a C4 transition between a state with the C-terminal helix and a state with the C-terminal loop, while the other modes showed a C2 transition between rotated states with two monomers with the C-terminal helix and 2 with the C-terminal loop. To extract C4 particles with the C-terminal helix: Set 1 particles with Mode 1 < −25, but Mode 0 and 2 between −50 and 50, Set 2 particles with Mode 2 < −25, but Mode 0 and 2 between −50 and 50. To extract C4 particles with the CTD loop: Set 1 particles with Mode 1 > 25, but Mode 0 and 2 between −50 and 50, Set 2 particles with Mode 2 > 25, but Mode 0 and 1 between −50 and 50. After combining the equivalent states from the two sets, there were 486,451 symmetry expanded particles in the state with the C-terminal helix and 486,812 symmetry expanded particles in the state with the C-terminal loop. 3D variability was run on each set again using the same parameters to further clean up the particles. For the particles with the C-terminal loop, the extremes of different modes all showed variation away from the average C4 state, so only the extreme particles (>50 or <−50) in each state were excluded. For the particles with the C-terminal helix, Mode 0 showed a transition between C4 particles with an intact pore domain and particles with a disrupted pore domain, while the other modes described a break from C4 symmetry. Only particles at the most intact end of Mode 0 were taken (>50), and the extreme particles from the other modes were excluded (> 50 or <−50). After reducing symmetry expansion, this resulted in 92,948 particles in the class with the C-terminal loop and 8,258 particles in the class with the C-terminal helix. These particles were then subjected to global and local CTF refinement and Nonuniform refinement with C4 symmetry. Finally, the particles were locally refined with C4 symmetry and map excluding the nanodisc density. The final pixel size was adjusted to 1.07 Å/pix based on comparison with a TRPV2 ARD crystal structure (PDB 2ETA) 55 . There were not enough particles in a stable state along the other modes of movement to be able to extract a good quality structure, and it is unclear whether these C2 and asymmetric transitions are functionally relevant.
For the 2-APB + CBD dataset, all data processing was done in cryoSPARC v 3.2. The 10,291 movies of the dataset were motion corrected using patch motion correction, binning the particles to 1.06 Å/pix. CTF fit of the motion corrected micrographs was estimated using patch CTF. This set of micrographs was then curated to remove suboptimal data, leaving 9354 micrographs. From this set, template-based auto-picking yielded an initial set of 1,709,112 particles. These particles were binned by 4 during extraction and subjected to two rounds of 2D classification to remove false positives and bad particles, reducing the particle count to 820,578. The particles were reextracted binned by 2 and subjected to two rounds of Ab-Initio Reconstruction, each time searching for four classes. The first round yielded 1 class with good TRP channel features, composed of 379,054 particles. The second round yielded two classes with good TRP channel features, composed of a combined 256,203 particles. These particles were reextracted to their final box size and refined by Non-uniform refinement with C4 symmetry and CTF refinement before symmetry expansion. These particles were then subjected to 3D variability analysis 27 , using a mask that excluded nanodisc density and searching for three modes of movement. The resulting modes were primarily assessed based on the cytoplasmic regulatory switch, as it is a well resolved region with a binary variable. Mode 0 represented a C4 transition between a state with the C-terminal helix and a state with the C-terminal loop, Mode 1 showed a C2 transition between rotated states with two opposite monomers with the C-terminal helix and two with the C-terminal loop, while Mode 2 showed an asymmetric transition between two adjacent monomers with the C-terminal helix and two with the C-terminal loop. To extract the two C4 classes, the 3D variability display in cluster mode searching for seven classes was used to analyze Modes 0 and 1. Two of these clusters corresponded to the most extreme ends of Mode 0, while excluding the most extreme ends of Mode 1. After reducing the symmetry expansion, this yielded 30,563 particles in the state with the C-terminal helix (the activated state) and 33,404 particles in the state with the C-terminal loop (the inactivated state). These particles were then subjected to global and local CTF refinement and Non-uniform refinement with C4 symmetry. The final pixel size was adjusted to 1.07 Å/pix based on comparison with a TRPV2 ARD crystal structure (PDB 2ETA) 55 .
Model building. For TRPV2 2APB_AC , the model of TRPV2 APO_2 (PDB 6U86) was used as a starting model and docked into the map in Coot 56 . For TRPV2 2APB_IAC , the model of TRPV2 CBD_2 (PDB 6U88) was used as a starting model and docked into the map in Coot. For TRPV2 2APB_CBD_AC , the model of TRPV2 CBD1 (PDB 6U84) was used as a starting model and docked into the map in Coot 56 . For TRPV2 2APB_CBD_IAC , the model of TRPV2 2APB_IAC was used as a starting model and docked into the map in Coot 56 . These models were then iteratively manually adjusted to the map and refined using phenix.real_space_refine from the PHENIX software package 57 .
The ligand restraint files for 2-APB and CBD were generated using the eLBOW tool from the PHENIX software package 58 . Due to the variable quality of drug density in the maps, 2-APB and CBD were only built into TRPV2 2APB_IAC and TRPV2 2APB_CBD_AC .
Pore profiles were generated using Hole 59 . Images of the models and maps for figures were generated using Pymol (Schrödinger), Chimera 60 , and ChimeraX 61 .
Whole-cell voltage clamp was performed with an EPC10 USB HEKA amplifier (HEKA Electronik, Lamprecht, Germany), much as previously described 31 . Signals were low passed at 1 kHz and sampled at 2-10 kHz. Patch pipettes were pulled from borosilicate glass tubes (TW150F-3; World Precision Instruments, Berlin, Germany) to give a resistance of about 3 MΩ. Cells were held at −60 mV and all recordings were performed at room temperature. The external solution contained: 140 mM NaCl, 5 mM KCl, 2 mM MgCl 2 , 5 mM EGTA, 10 mM glucose, and 10 mM HEPES, pH 7.4 (adjusted with NaOH). Calcium was omitted in order to avoid desensitization. The standard pipette solution contained: 140 mM KCl, 2 mM MgCl 2 , 5 mM EGTA and 10 mM HEPES, pH 7.4 (adjusted with KOH). A gravity-driven glass multi-barrel perfusion system was used to bath apply solutions. Data acquisition and off-line analyses required Patchmaster/Fitmaster software (HEKA Electronik, Lamprecht, Germany) and Origin 8.5.1 (Origin Lab, Northampton, MA, USA).
Rosetta docking of 2-APB. The starting structure of TRPV2 2APB_IAC was refined using Rosetta relax with electron density constraints (see Rosetta scripts and command lines used in Supplementary Note 1). The initial fitted 2-APB molecule was removed for the refinement. We generated 1000 decoys and selected the topscoring model for the docking step. OpenEye OMEGA (OpenEye Scientific Software) 62 was used to generate conformers for 2-APB molecule. The primary amine group on 2-APB was protonated to be consistent with its form at the physiological pH. We used RosettaLigand for docking of 2-APB to the top refined model of TRPV2 (see Rosetta scripts and command lines used in Supplementary Note 2). The details of the RosettaLigand docking algorithm have been previously described 29,30 . Here, we used a sampling radius of 3 Å to perturb 2-APB around the binding pocket identified by the electron density. An initial random perturbation of 0.5 Å was added before the Monte Carlo docking cycles to further randomize the docking process. A total of 5000 docking models were generated. Rosetta ligand binding energy represented by the interface_delta_X term was used to select the top docking model. All-atom RMSD was used to compare all docking models to the initial fitted model ( Supplementary Fig. 5).
Reporting summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The atomic coordinates and cryo-EM density maps generated in this study have been